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ABSTRACT 











Load handling technology in the sea environment to place heavy and bulky 
ads accurately on the ocean floor is in its preliminary stages of development. 
e load handling system is usually mounted on a platform which interacts with 
e sea surface, imparting dynamic loads to the lifting lines that are not ex- 
rienced on land. This study evaluates the response of a specific load handling 
rystem composed of a crane and constant-tension winch mounted on AMMI pontoon 
arge to sea states 3, 4, and 5. The motions of heave, pitch, and roll were 
oupled by phasing at specific load suspension points and the rms values of the 
erticle displacement and accelleration at these points were calculated using 
pectral analysis techniques. It was found that the constant tension winch could 
mpensate for the motions of sea state 3 but that the probability of exceeding the 
inch capability in sea state 4 was patee enough to limit operations to missions of 
ery short duration. 
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1. INTRODUCTION 


Construction involves the lifting and transportation of heavy equipment and 
all kinds of materials. For dry land application, a whole spectrum of simple and 
complex load handling systems have been designed and built to achieve the comple- 
tion of construction projects of all kinds and magnitudes. As man now moves to a 
completely different environment for construction, the ocean floor and the sea 
above it, it is found that load handling technology is in its preliminary stages of 
development at best. To accomplish underwater construction there is a need to 


lift and lower and also recover heavy and bulky objects in the open seas. 


One of the conventional systems often used for this type of load handling 
consists of a winch and a line with the line being run out on some kind of boom. 
The winch and boom are mounted on a small floating platform and a load is sus- 
pended at the end of the line. As the platform interacts with the waves of the ocean 
surface, the platform motions cause the load to oscillate and this in turn generates 
oscillation in the line tension. In order to design a load handling system for a 
specific job, the ocean engineer must be able to predict how the dynamic tension 
in a line from a given platf{« rm will vary with the sea state and specific loads at- 


tached to the end of the lin::. 


Because of the U.S. Navy's interest in lifting and lowering loads to 6, 000 
foot depths, a development program has been carried out at the Naval Civil Engi- 
neering Laboratory (NCEL) to predict the dynamic stress response of lifting lines 
in Oceanic operations. These ainhes eo and others have resulted in a method 
of predicting lifting line stresses for various ideally shaped loads using the sea 
induced motions of the platform at the last point of contact with the line as one of 


the major inputs. 


Much work has been accomplished in hydrodynamics for the prediction of 
ship motions in both regular and irregular seaways since the introduction of the 
strip theory of ship motions by Korvin-Kroukovsky in 1955. ‘ Rightfully this re- 
Search into the optimization of the hydrodynamic characteristics of transportation 


type ship structures belongs in the field of Hydrodynamics or Naval Architecture. 








10 
The Ocean Engineer on the other hand is most likely concerned with a working 
platform that is barge shaped and may be moored or remains fairly stationary in 


the seaway. 


Thus, this thesis will concern itself with the prediction of the motions of 
a fairly simple ocean engineering construction load handling system in a seaway. 
The goal of this work is to provide a simplified method of analysis that will give 
a good first estimate of platform vertical accellerations at points on the platform 
where lifting lines would receive the greatest effect from the barge motions. The 
data developed by this study should be accurate enough to allow sizing of lifting 
lines within engineering safety and prediction of the operational limitation of the 


system as a function of sea state. 
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2. SCOPE OF THE THESIS 


In 1970, the Navy executed project AFAR? which involved the implantment 
of an acoustical array on a sea mount in 1300 feet of water. It was determined 
that the most Economical load handling system to support the project was to con- 
struct work barges from 90 feet x 28 feet x 5 feet AMMI pontoon sections or from 
basic 5 feet x 5 feet x 7 feet pontoon cans. In the operations alan the construc- 
tion was administratively limited to sea state 3 on the Beaufort scale, however no 
analysis was provided to indicate what magnitude of barge motions or line tensions 


that could be expected at that or higher sea states. 


This thesis will specifically analyze the interaction of the Crane Barge of 
Project AFAR shown in Figures 1 and 2 with various sea states starting with sea 
state three. The barge is composed of one AMMI pontoon section fitted with a 
Harbor Master, 180 HP, propulsion unit in each corner, a 30 ton capacity stand- 
ard construction crane, and a Leithiser constant-tension winch. The propulsion 
system can be controlled so that the barge can remain on a specific station, with- 
out mooring, within the accuracy required for the implantment. The physical 


characteristics of the barge are listed in Table 1. 


The method of analysis to be used will employ the linearized equations of 
motion for pitch, heave and roll for a regular unidirectional train of waves ap- 
proaching the barge from the starboard quarter. The wave frequency dependent 
coefficients and the wave excitation amplitudes and phase lags of the equations 
will be determined by comparing a method of calculation for the coefficients pub- 
lished by Muga of NCEL to the experimentally determined coefficients of Vugts 
for a rectangular cylinder with a beam~-draft ratio of 8, and using the coefficients 
felt to be most accurate. The equations will then be solved for input waves of 
seven different wave lengths. A response operator for each of these wave lengths 
for the verticle acceleratim will be determined with phase coupling of pitch, 
heave and roll at four points of interest: the deck amidships at both the bow and 
stern and the end of the crane boom with the boom at both the port beam position 


and 60 degrees after the port bow. A smooth curve for the response operator 
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Figure 1. Project AFAR Crane Barge 
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Figure 2, Drawing of Project AFAR Crane Barge 
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PHYSICAL CHARACTERISTICS OF AMMI 


Length 

Beam 

Calm water draft 
Displacement 
Weight 


GM 


Mass 


CRANE BARGE 


90 feet 

28 feet 

3.0 feet 

202 L. Tons 

564, 480 lbs 
147.25 feet 
(owes 10° Slug-ft" 


7.125 x 10° Slug-ft 


17, 640 Slugs 
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at each of the points will be drawn and using spectral analysis techniques, a re- 
sponse spectrum for several sea states will be determined. From the response 
spectrum, the R.M.S. value of the accelerations for each sea state will be calcu- 
lated and statistical information developed which can provide the basis for limiting 


operations above certain sea states. 





—_ 
- 
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3. METHOD OF ANALYSIS 


3.1 Equations of Motion 


The linearized equations of motion for the barge will be developed along 
the technique suggested by ROMO. The motions are separated into motions in 
the horizontal plane and motions in the verticle plane. Because of the difficulty 
and complexity in determining cross coupling coefficients and the uncertainty in 
the state of the art in handling all six degrees of freedom, only pitch and heave 
will be considered in the verticle plane and only roll in the horizontal plane. The 
excitation due to waves is the only dynamical response term that will be considered 
on the right hand side of the equations. It is assumed that the barge can maintain 
position with no forward speed and thus control and propulsion forces are neglected. 
The regular wave train is assumed to be unidirectional and approaches the barge 
from the starboard quarter or 45° off the bow. To simplify the analysis, the reg- 
ular wave is broken down into two equal components of wave length V2y, one from 
directly head-on and the other from beam-on. Pitch and heave will be affected by 
head-on waves while roll and heave will be affected by beam-on waves. The cen- 


ter of gravity is assumed to be at the water line, amidships. 


The notation to be used is that D = < and a symbol with a dot over it, 6, 
is the first derivative with respect to time. A symbol with two dots over it, z is 
the second derivative with respect to time. The derivatives for the hydrodynamic 
oe) will be written as ae which is the partial derivative with re- 
spect to the vertical accelleration taken at the original conditions with the barge 


coefficients, ( 


at rest. In this notation it should be remembered that p = ¢ and q = 6. 
Thus the general linearized equations of motion for the barge system are: 


Heave 
2 2 1(wt -€7) 
_ + + . + mx + + = - 
| (Zs mv) 1D) eZee) Z.,_| Z (2, mx.) D ae Z ; ] Ss) Ze 


fa] n+ Cay] O= At) 


a 





-.t- 








ut) 
Pitch 


2 va 
|(M, + mx )D° +M D+M Jz +|@, -I)D-+MD+M |@=-M oilwt~ ey) 
Ww G W Z q gy q 3) il 


Roll 


_ 2 _ i(wt -€ 1) 
Lk, L)D'+K D+ K, 2 =-Ke (3) 


where Z 1? M, and K are the amplitudes of the excitation and ¢ is the phase lag. 
Since there is no Cee speed, w was UE in place of W , the frequency of en- 
counter. Thus Weta ow the wave frequency. The eee on the left 
hand side of the ee, are independent of time but may be a function of w. Due 
to the choice of the CG and the constant section of the barge, the coupling coeffi- 


cients will be zero. 


3.2 Solution of the Equations 


Solution of the equations of motion is accomplished by solving for the roots 
of the determinant of the coefficient matrix for equations (1) and (2) or by just 


solving for the roots as in equation (3). The solution will be of the form: 


mm 
_ Ot i(wt - €;) 
f(t) = ) Oe +f), @ j (4) 


§ where T(t) is the chosen independent variable, m is the number of independent 


' variables in the system, C, arc the constants of integration, o, are the roots of 


k k 
the determinant, ues is the amplitude of the steady state response and ¢ olate 


phase lag. 





flaw 
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The solution of equations (1) and (2) in determinant form is: 


A,(t) ay, a, At) 
A(t) a a. A(t) 
2 22 21 “2 
LO) ara Oa 220); - = = (6) 
a a a a 
12 mae 
451 90 Son = 498 


The steady state solution of equation (8) is: 


K, ol (wt -eK - 1) 
Ra as (7) 
Ke (ee yar Pe eK wy" 
Ky p x p 


where 


l K WwW 
T = tan ae (8) 
3 = PS ence 


| 3.3 Evaluation of the Coeffieients 


The hydrodynamie eoeffieients are independent of time beeause they are 
made up of various derivatives defined by referenee only to the initial equilibrium 
condition. They may depend on the vehiele size, shape and inertial distribution 
and on the frequeney of exeitation sinee the motion of the vehiele on the free sur- 
face generates waves. The frequeney of this response must be the same as that 
| of the exeitation in order for the eoeffieients to remain time independent. A sys- 
tem of axes is defined for this analysis as shown in figure 3 with the forees X, Y, 
and Z positive in the direetion shown, piteh positive in the bow up direetion and 
roll positive with the starbord side down. The origin is taken at the center of 
gravity of the barge whieh is assumed to be direetly above the keel amidships at 


the ealm water line. Thus x G = CG ae, G = 0. Eaeh of the eoeffieients will 
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generally include derivatives which are cither inertial, damping, or hydrostatic 
in nature. In the case of damping, only damping duc to wave generation will be 


considered. 


There has been significant research carried out in recent years to develop 
methods of theorctical calculation for the more important hydrodynamic coeffi- 
cients. References 7 through 12 are representative of the current state of the art 
and provide an excellent reference list. Coefficients for motion in the verticle 
plane are better defined and understood than those in the horizontal plane. Some 
of the coupling coefficients have never been measured due to the difficulty in de- 


vising an accurate experiment. 


Methods for calculating the coefficients were developed for the Navy by 
Kaplan and Putz (1962)"° for the CUSS I Project MOHOLE Barge in a moored 
condition and by Muga (1966)° for the missile recovery barge FISHOOK, also in 
a moored condition. Although the barges considered in these studies were slightly 
more complex, the method of coefficient calculation is applicable. Frank (i967)** 
has published section: coefficient values for the center section of a series-60 
ship, however, with a smaller beam-draft ratio than the barge of this study. 

Vugts (1968) ‘ has also published the sectional coefficients of long rectangular 
cylinders for various beam-draft ratios, some closely approximating the condi- 


tions of this study. 


To simplify this analysis, the method of calculations published by Muga® 
and the experimental data of mee. will be used as best fitting a rectangular box 
shaped barge. The sectional derivatives calculated by Muga's equations will be 
compared with the experimental data of Vugts and the experimental derivative will 
be used whenever available. Because of the symmetry of the barge, it is assumed 
that the sectional coefficients are constant over the length of the barge. Details 
of the derivation of each coefficient is presented in Appendix A. Table 2 lists the 
nondimensional values of the frequency dependent coefficients calculated by Muga's 


equations compared to those determined hy Vugts. 
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TABLE 2 


HYDRODYNAMIC COEFFICIENTS 


(Calculated vs Experimental Data) 
(Note: p = sectional mass density) 





d 50 100 «=—ss150,— 200 300 5.00 600 
a 
Heave Added Mass =—~ 
pA 
fiea’ Zo 2.76 Det 4,27 o. 01 een 5.41 
Vugts’ 2.00 2h US 3.08 Be 0 3.90 4,50 5. 00 
Z 
| Heave Damping = — P 
pA 2g 
Muga 0.51 2.30 3.24 3.61 3.92 3.09 2.85 
Vugts 1.45 2.13 2228 VA Pag 1h AY) 2.10 1.95 
M. 
Pitch Added Moment = —2 
pAL? 
Muga 0.209 0.230 0.312 0.356 0.418 0.439 0.451 
| Vugts 0.229 0.229 0.256 0.308 0.325 0.375 0.417 
Pitch Damping = mar ie 
pAL 5 
Muga 0.043 0.191 0.270 0.301 0.293 0.257 0.237 
Vugts Do Leah 0.177 0.187 0.185 0.183 0.175 0.162 
K. 
Roll Added Mass = —P- 
2 
pAB 
Muga (insufficient information to make the calculations) 
Vugts 0.098 O22 0.139 0.146 0.148 0.146 0.145 
» [B 2 
Roll Damping = Pp j= (Note: Muga Values all x 10 “) 
2 j2¢ 
pAB 
Muga 0.0153 0.0104 0.0081 0.0065 0.0050 0. 0029 0.0023 


Vugts 0.0300 0.0375 0.0270 0.0180 0. 0080 0. 0050 0.0030 








21 


Figure 3. System of Axes 





t 


TABLE 3 


HYDROSTATIC COEFFICIENTS 


Z. - Mieise0> Int/it 
8 

Mo = 1,086 x 10° Ibf - ft 
7 

K ~ 8.34x 10° lbf - ft 
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3.4 Wave Excitation 


The wave induced effects on the barge are composed of inertial, damping 
and displacement components of the wave induced force acting on individual ship 
sections of length. For this study it has been assumed that a unidirectional train 
of regular waves approaches the barge at a 45 degree angle with the bow at some 
wave length, }. The incoming wave length has been broken down into a head-on 
component and a peameon component, each of wave length, ./2\. The wave excita- 
tion forces and moments on the right hand side of equations (1), (2) and (3) will be 


defined in terms of the hydrodynamic coefficients previously derived in Section 3.3. 

The details of the definition of the wave excitation forces and moments are 
explained in Appendix B. It should be noted that the general form of the equations 
for the forces, an or moments, M in Appendix B are of the general form: 


F or M =A cosut+C_ sin wt (9) 
W W W W 


which can also be written as 


if ene 1 el : +C : COS (Wt - €) (10) 
Ww W W Ww 
where % oy 
= tan — 11 
€ =tan A (11) 
Ww 


If it is understood that only the real part of complex notation will be considered, 
then equation (10) can be written as 


Home Ac oe te °) (12) 
W W W 


W 
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Comparing this expression to the right hand side terms of the equation of motion, 


it can be seen that for heave as an example. 


ee = A af ae 
Z Oe (13) 


and 


ai(ut ez) _ i(ut- ¢) 14) 


Evaluation of the excitation forces in this form gives the maximum ampli- 
tude of the excitation for one wave and the phase relationship of this amplitude com- 
pared to the wave with its crest exactly admidships at the center of gravity. The 
excitation amplitudes Z 


1 1 1 © 
calculated and are tabulated in Table 4 for each of the frequencies of interest. For 


, M, and K, and the phase angles oo € ande 5 have been 


comparison, the values of maximum generated wave amplitudes and phase angles 
7 ae 
measu:s! by Vugts as a consequence of a known excitation of his model are also 


tabulat-:i in Table 4. 


3.9 Response Operators 


The equations of motion can now be Solved for the desired variables 
Z ) = 45 degrees, © and ¢ for any wave length or frequency of regular wave. 
(Note: w is the angle between the bow of the barge and the incoming component 


of the wave train.) The solution vectors will be of the form; 


nee R. cos (wt + oF (15) 


or 


ge ty dr) 


12 (5) R. (16) 


where r represents any of the motions, R. is the magnitude of that motion or am- 
plitude response operator, and .. is the phasing of the maximum possible motion 
relative to maximum positive wave displacement at amidships. O.. is sometimes 


called the phase response operator. Thus for a unit wave amplitude, r(t) is the 
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response operator which is frequency dependent. The values of the motion re- 
sponse operators for heave, pitch and roll are listed in Table 5 and plotted as a 


function of w in Figures 4, 5 and 6. 


Having calculated the motions, the vertical displacement, Z_, velocity 


P’ 
) can be calculated by 


: a i 
Zp: and acceleration, Zp p’ Yp 
coupling the motions from each component of the wave. Thus the vertical re- 


of any point P = (x 


Sponse operator of any point can be written as: 


Zp =| 4 450 Op) + Op) | (17) 
tA = 1Kh [, 45° 7 OK) + op) | (18) 
Z, =ow 4, 45° ©%p) +6 vp) | oe 


It must be remembered that each of the response operators is a vector and thus 
must be added keeping track of all the phase response operators 6. Figure 7 is 


-an example phase diagram of this type of vector addition. 


Since two of the four points of interest in de.vrmining maximum accelera- 
tion includes a boom, the boom can be included in the response operators of equa- 


tions (17), (18) and (19) by merely letting x_ =1 cos acos y and Yp= 1] cos @ Sin y 


Pp 
where | is the boom length, a is the elevation angle and y is the azimuth angle. 
Values for the response operator at the four points of interest are tabulated in 


Table 6. 
3.6 Description of the Seaway 


Since the response to any given wave frequency Can now be calculated, it 
is necessary to apply the frequency spectra of an irregular sea using the energy 
density technique described by Loukakis (1970)! to predict the statistical proper- 
ties of ship responses in a seaway. The spectral family that will be used to 
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1 
describe the seaway is the Pierson and Moskowitz (1963) ° family for a fully 


developed seaway. The general form of the spectral density function is given by 
-5 —4 
Q(w) =aw exp(-Bw ) (20) 


2 
where Q(w) is the one sided spectral density function with units (ff x sec). The 


parameters a and § for a fully developed sea are given by: 


Z 
a = 0.0081 ¢g Cryoorl 
(21) 


8 = 0. 0324 oa jay” cel | 


a 
where H® is the significant wave height which by definition equals four times the 

root mean square (rms) of the spectrum. Specific values for the spectral density 
function for three sea states at the frequencies of interest are tabulated in Table 7 


along with corresponding wind speeds. 
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3.7 Response Spectra 


All absolute and relative barge motions, velocities and accelerations are 
nearly linear with respect to the wave. Thus the barge behaves approximately as 
a linear, time invariant filter and the output of the filter will have the same distri- 
bution as the input. If the sea spectral density function is plotted as proportional 
to is where a is the wave amplitude and 6w is an incremental frequency element 
and the barge response operator is plotted as proportional to (z/ay-, both as func- 
tions of w, it is clear that by the superposition principle: 


2 2 


Z — 
x = (22) 


|” 
pe) 
[bo 


When ra /5 is plotted as a function of w, a response spectrum results. A quantity 
is now defined which is the average value of the response average over the fre- 
quency range. This average is obtained by integration over frequency or is the 
area of the spectrum. The root mean square, rms, of the response amplitude is 
then, rms =,/area, which represents a statistical average of the irregular random 
responses. From the rms value, other statistical averages such as the average 


value of the 1/3 or 1/10 highest response can be calculated by the relationship: 


Response ’/N = £(N) (rms) (23) 


The values of N and t(N) which are of interest in this study are listed in Table 8. 


The values for calculated responses at the points of interest for the seven 
wave lengths investigated are tabulated in Table 9. These values were plotted and 
the area under the curve determined graphically. The rms values for each of the 


responses are listed in Table 10. 
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TABLE 4 


EXCMIAlION AMPLITUDES AND PHASE ANGLES 


Ww Lely = wea Z elGre 00 WeoGeosos  O200 


3 
Amplitudes (Forces - lbf x 10°, Moments lbf-ft x 10 ) 


Z 7.84 34.5 [eel 1833 108.5 128.8 130.3 
e= 45~ 

Zz (Vugt) $1.5 82.1 96.6 107.8 116.0 124.0 128.7 
M. 291 1610 VAC Re RCRsT:: iZodm724 meno24 
Ky 96 265 287 255 196 125 106 
K (Vugts) 200 319 302 266 DOG. 128 109 


Phase Angles (Degrees) 


ez, # = 45° Ins Beye 216 212 ce oS) “193 
€, (Vugts) 135 is 40 30 Zo 20 10 
eu 24 S22 308 302 293 235 283 
f 90 270 SG 270 270 270 270 
€ (Vugts) 0 270 225 (remainder not measured) 


K 





TABLE 5 
MOTION RESPONSE OPERATORS 


a) 220) W420 “1216 00M 0.82 


Heave (ft /ft) 


z/a 0.063 a7 07 0 0.840 0.943 


Pitch (Radians/ft x 10> 


O/a 4.36 18.42 16.75 Wa ly 


Roll (Radians/ft x 1079 


P/a 1.28 2A BG 3.16 2.40 


OF 


OSS 9 


i 


lis 


63 


ie 


52 


28 


OSS 


OF 999 


1.29 


z/a (f{t/ft) 


(Radians/ it} 


©/a x 10 


Figure 4. Heave Response Operator 
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Figure 6. Roll Response Operator 
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Figure 7. Sample Phase Diagram 
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TABLE 6 


for ONoE OPERATOR Al POINTS OF INTEREST 


G) fp NI 1 £42 1.16 1.00 
Bow Amidships (x =45 ft) 
EE ° eee 
e/a (ft ft) 0) Aw 0.92 0.63 0.71 


n. 2 
z/a(fit/sec -ft) 0.90 1.84 Osis Oya 


Stern Amidships (x==45 ft) 
TL EEE 


z/a (it/ft) 0.16 eZ ls OF 130 


bya (ft/sec’-f)0.65 2.05 1.84 1.30 


Boom at Port Beam (Boom = 60 ft,a=60°, ¥=-90°) 


z/a (ft/ft) Oa, 0.40 0.64 0.74 


Bya(it/sec’-ft)0.42 0.81 0.86 £0.74 


Boom Aft Port Bow (Boom = 60 ft, X= 60°, ¥'=-60° 


e/a (ft/ft) 0.10 0.23 0.54 0.66 


Ma (ft/sec--fi)0.41 0.47 0.73 0.66 


OZ 


0.89 


0.60 


.03 
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Aa), 


~13 


~45 


or 


Ae, 


cele: 
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Sl 


Bees. 


Bly 


AS): 


709 


oe 


Ache: 
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Boe. 
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TABLE 7 


PIERSON-MOSKOWITZ SPECTRAL DATA 
* FOR A FULLY DEVELOPED SEA 


Sea State Data 


Sea State 3 4 5 
Wind Speed (KTS) 14.5 19.5 23.2 
H® (ft) 4 7 10 

rt | 8.294 8.294 8.294 
B 2.07 0.677 0.332 


Spectral Density Function (ft”-sec) 











w Q,(w) Q, (w) Q, (w) 
0.58 0.0018 O2323 6.77 
0.634 0.0191 I AY 10.50 
0.82 0.228 5.02 10.75 
1.002 1.06 4.18 5.88 
1.16 1226 DA Pas 3.30 
1.42 0.865 are 1a) 
2.01 O72726 0.242 0.247 
TABLE 8 
1/N 
VALUES FOR: RESPONSE =f, (N) (rms) 

N £,() 

1 1,25 

3 2.00 

10, 2.55 

10 Oeoe 


10° 5.27 
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TABLE 9 

RESPONSE SPECTRUM VALUES 
fo 2.01 1.42 ello 1.00 0.82 ORiG3 OPES: 
sea State 
Bow Amidships - Displacement 2” /Sw Go 
3 e011 On 73 0.50 Ues5 02-210 0.018 OG 
4 Orso 2 02 ON, Zeeks 4.48 souls! Dod 


° OG elke le el eS) G00 Dey 6.37 


Bow Amidships - Acceleration ve /Sus (ft /sec”) 
3 O) Saks 2.94 Deol 0.595 0.08 0.003 0.000 
4 O) calle 4.19 1.90 Zales F180 Oe em 0.035 


° O20 4.52 Dra, IS Sas ee, 0.724 


Stern Amidships - Displacement 2” /su 


= 0.006 0.89 1) By, 1.80 0.37 0.024 0.002 
™ 0.006 1.26 5.12 7 ao 8.06 oz 0.39 


. 0.006 es? O20 a 29)7, 17.30 oe 8.10 


otern Amidships - Acceleration 2% /Su0 
mm 0.090 3.62 4.26 tee rele 0.004 0.000 
m™ 0.101 Se 10 SZ TER ARE 3.64 0.258 0.044 


5 0.103 5.56 11.16 oo Te 9 Z2UGIS O29 16 
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TABLE, Ge(Cont a.) 
RESPONSE SPECTRUM VALUES 

mo 2.01 1.42 ier s' O10 Orse2 O263 0.58 
Boom at Port Beam - Displacement Bo Ned 
mB 0.092 0.14 #O.51 40.55 0.18 0.018 0.002 
4 0.003 0-40 aa Desh 3.94 ees 0.311 
9 Oe 03 O22 1.34 $208 8.45 sei G02 
bea State 
Boom at Port Beam —- Mecelérstion’:?/$w 
3 0.040 O57 O92 0655 O05] 0.002 0.000 
4 0.043 Ute ZOO Ag aes: Oil 9 O2035 


9 0.044 0.88 2.42 3.06 ee leuk L512 0.74 


Boom Aft Port Bow - Displacement 2° /S$u9 


3 0.002 0.046 0.37. 40.47. #20.15 0.017 0.002 
™) 6.003 0.066 0.79 41.84 &2«1|3.38 #+#+.1.09 ~° &0.29 


o 0.003 0.071 eo / A Bese PAs oe hes B07 


Boom Aft Port Bow - Acceleration 22/5 
3 O03 7 Os, Orb” 0.47 0.069 0.003 0.000 
4 0.040 E20 1.43 1.84 Re: O76 O30 3 


9 0.041 OE Zs) 1.74 2ao6 Soy 1,49 Oates: 





Location 


Bow Amidships 


Stern Amidships 


Boom at Port Beam 


Boom Aft Port Bow 


TABLE 10 


RESPONSE RMS VALUES 


Motion 


72, Vile 


7, ft/sec 


loo 


UW gs) 


24 


1.09 
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0.53 


On71 


0.44 


Or / 


es 
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pon 


.04 


OW 


mZO 
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4. Discussion of Results 
4.1 Response Operators 

The response operator amplitudes for heave, pitch and roll have 
been presented in Table 5 and Figures 4, 5 and 6. The heave operator 
behaves as expected approaching a maximum value of the wave ampli- 
tude. The pitch operator is slightly lower than the maximum wave 
Slope at the low frequencies. If the maximum wave slope is given 
by, then at \= 600 ft, the slope would be only about 1.05x10 we 
(Radians/Unit amplitude). The values of the operator are in general 
agreement with the values found by Kaplan and Bite for CUSS I. 
The roll operator behaves as ecpected but is about an order of 
magnitude lower than values determined by Kaplan and Putz. It 
should be mentioned that CUSS I has a more ship type hull than the 


7. 
crane barge. The roll coefficients of Vugts also had considerable 


spread in the data which may account for the difference. 


4.2 Response Operators at Points of Interest 

The response operator amplitudes for vertical displacement 
and acceleration for the four points of interest have been presented 
in Table 6. The values for bow and stem amidships reflect the 


coupling of pitch and heave. The contribution of roll to the motion 





o7 
of the crane boom was fairly small compared to pitch and heave 


due again to the small value of the response operator. 


4.3 RMS Values 
The results of the spectral analysis were presented in Table 10 

in the form of rms values for vertical displacements and accelerations. 
The value for displacement in sea states 3 and 4 is felt to be accurate, 
however sea state 5 had significant contribution between 1.42 and 2.01 
radians where the shape of the curve was not calculated. The acceler- 
ations are more accurate however since the values at both ends of the 
frequency span chosen are negligible. All of the rms values would have 
been more accurate with more eee especially at the low end of the 


frequency scale and betweenW= 2.01 and W= 1.42. 


4.4 Limitation on Operations 

Now that the rms amplitude values of verticle displacement and 
acceleration are known, the probability of having a total response of 
any magnitude can be calculated from equation (23). For example, for 
sea state 4 the rms of the displacement amplitude at the bow is 1.35 
ft. Using Table 8, for 66.7% of the time, the displacement would be 
less than about 2.7 ft. Only one displacement in 100 would exceed 
4.5 ft and one in one million would exceed 7.1 ft. Since the average 


period of sea state 4 in this study is about 5 seconds, this means that 
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the occurrence of an excitation which would produce a displacement response 
greater than 7.1 ft would be only once ina time period of about 58 days. 
Based on this type of analysis, the displacements and accelerations re- 
ceived as a result of sea state 4 are too large for anything but very short 
mission duration weight handling. Sea state 3 rms values of this study 
should be satisfactory for operations when synthetic line and the leithiser 
constant tension winch are used to compensate for part of the displace- 


ment and acceleration. 


5. Conclusions and Recommendations 

It appears that the spareeonell limitation on Project AFAR to sea 
state 3 on the Beaufort scale was a bit conservative since this corresponds 
to a wind speed of 10 knots whereas the Pierson-Moskowitz sea state 3 of 
this study corresponds to a wind speed of 14.5 knots. Operations in this 
latter sea state should be acceptable if the Leithiser constant-tension 
winch is used with synthetic line to compensate for the motions. For sea 
state 4 of this study, the probability of exceeding the capability of the 
weight handling system is great enough that operations would be limited 
to very short load lifting missions at best. 

It is recommended that this method of analysis be reduced to a 
computer program for all heading angles at more frequency points to pro- 
vide more accurate rms values. It is also recommended that more research 
be done on refining the hydrodynamic coefficients for roll for shallow draft 


box shaped hulls to eliminate the uncertainties in the data of Vugts. 





nO. 


ol. 


12, 


13. 


14, 


39 
REFERENCES 


. NCEL Report TR-543 "Mechanics of Raising and Lowering Heavy Loads in the 


Ocean: Experimental Results", Muga, B.J. Sept 1967. 


. NCEL Report TR-652 ''Feasibility Study and Comparative Analysis of Deep 


Ocean Load Handling Systems" by Davis, D.A. and Wolfe, M.J. Dec 1969. 


. NCEL Report TR-703 "Dynamic Stress Response of Lifting Lines for Oceanic 


Operations" by Liu, C.L. Nov 1970. 


. Korvin — Kroukovsky, B.V. "Investigation of Ship Motions in Regular Waves," 


Transactions of the Society of Naval Architects and Marine Engineers, Vol 63, 
pp. 386-435 (1955). 


. Naval Facilities Engineering Command, Project AFAR(AUTEC), Management / 


Operation Plan 1/27/70. 


NCEL Technical Report R-440, ''Hydrodynamic Analysis of a Spread-Moored 
Platform in the Open Sea’, Muga, B.J. Aug 1966. 


. Laboratorium Voor Scheepsbouwkunde Report No. 194 'The Hydrodynamic Co- 


efficients for Swaying, Heaving and Rolling Cyclinders in a Free Surface" by 
Vugts, J.H. Jan 1°58. 


- Abkowitz, M.A. "Stability and Motion Control of Ocean Vehicles" MIT Press 


(1969). 


Abkowitz, M.A., Vossilopolous, L.A., and Sellers, F.H. "Recent Develop- 
ment in Seakeeping Research and its Applications!, Trans SNAME, Nov 1966. 


Mandel, P. "Principles of Naval Architecture — Ship Maneuvering and Control", 
SNAME 1967. 


Korvin — Kroukovsky, B.V. ''Theory of Seakeeping'" SNAME 1961. 


Lewis, E.V. ''Principles of Naval Architecture — The Motion of Ships in Waves", 
SNAME 1967. 


Kaplan, P and Putz, R.R., 'The Motions of a Moored Construction Type Barge 
in Irregular Waves and Their Influence on Construction Operations", NCEL 
Contract NBY-32206, August 1962. 


NSRDC Report 2375, "Oscillation of Cylinders in or Below the Free Surface of 
of Deep Fluids", Frank, W., October 1967. 





15. 


16. 


7. 


40 


REFERENCES (Cont'd. ) 


MIT Report 70-3, "Computer Aided Prediction of Seakeeping Performance in 
Ship Design", Loukakis, T.A., August 1970. 

Pierson, W.J. and Moskowitz, L., "A Proposed Spectral Form Based on the 
Simularity Theory of S.A. Kilaigorodski", Geophysical Sciences Laboratory, Re- 
port No 63-12, N.Y. University, (1963). 


Vossers, G. (1960) "Fundamentals of the Behavior of Ships in Waves", Inter- 
national Shipbuilding Progress, Vol 7, No. 65, Jan 1960. 





4] 
APPENDIX A 


EVALUATION OF THE HYDRODYNAMIC COEFFICIENTS 


Al Inertial Coefficients 


The coefficients of the accelleration terms are the result of the hydrody- 
namic forces and moments caused in reaction to the accelleration of a body through 


the water they are negative in value and thus add to the mass a inertia of the body. 
Al.1 Heave Added Mass, Zs and Ds 


Z. is the added mass in heave which is in phase with the verticle acceller- 
ation in equation (1) while as is the added inertia in phase with the angular accel- 
leration of pitch. Muga® defined these coefficients over a unit section of ship 


length as; 


2 
* 
Zz, =-(7> See ae 
Ww 2 4 


(24) 

bs 
mee ey tdt 95 
4 Js W (29) 


Ss 


where B* is the sectional beam, C is a frequency dependent correction factor, 
and € is a dummy variable along the x-axis which is integrated from the stem, 
Eg? to the bow a The limits of integration will be deleted in future integrals 
unless they differ from the above. It should be noted that due to Z being con- 
stant for each section, the coupled pitch-heave added mass, bee will be zero for 
all frequencies. The values calculated by equation @4) for Z; are tabulated in 


7 
Table 2 along with the comparable experimentally determined values of Vugts. 
Al.2 Pitch Added Inertia, M.. and ara 


These derivatives are in phase with the verticle accelleration and pitch 


: 6 
angular accelleration and have been defined by Muga as: 





42 


M. = \ Zz; €dé (26) 
M.= (Zz. é7dé (27) 
‘ee 


For the same reasons as expressed for a in equation (25), M.. is also equal to 
zero. The values of te calculated from equation (27) are listed in Table 2 along 


with the experimental values of Vugts. ! 
A1.3 Roll Added Inertia, se 


Muga calculated =a from the effective radius of gyration, KO by the rela- 


tionship: 


mk “=(I_ -K.) (28) 


The effective radius of gyration was calculated from the experimentally deter - 


mined roll natural period of FISH HOOK by the relationship: 


Drea (29) 
[g(GM)]* 


Since the natural roll period of the AFAR crane barge is not known with any accu 


racy, only the experimental values of a determined by Vugts are listed in Table 2. 


A2 Dynamic Damping Coefficients 


The hydrodynamic damping forces and moments are in phase with the body 
velocity and are proportional to the dissapated wave energy. Viscous effects of 
the fluid have been neglected in this study. Thus the damping coefficients are the 


coefficients of the velocity terms of the equations of motions. 
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A2.1 Heave Damping, Z and an 


Muga” derived the heave damping coefficients as follows: 


2— 2 
ps (A) 
2. ee \ perercuenm Se 
W 
2— 2 
ps (A,) 
a  edtic..=1 (31) 


W 


where A, is the ratio of the amplitude of the heave generated two dimensional waves 
to the heaving motion of the section and Cz and Cc. © are three dimensional damping 
factors. Again because of the constant sectional coefficient, a is equal to zero in 
equation (31). The values of Z calculated by equation 80) are compared with the 


experimental values of Vuzts ! in Table 2. 


A2.2 Pitch Damping, M and ei 


These cocfficients in phase with the verticle velocity and the pitch velocity 


are defined by niyo as: 


— 2 
ps(A,) 
a eens? oa : = = bye 
M coe 3 fdé;C, =C_ 67) (32) 
— 2 
== oo ee : — Spb. 
a “oe ) ie 2) SB eys ee ) 


S is once again equal to zero due to the coefficient being constant over the length 
{the body. The values of td calculated from equation (33) are compared to the 


xperimental values of Vugts in Table 2. 
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A2.3 Roll Damping, a 


Muga derives the roll damping coelficient as: 


_— a foil 2 2 
K=-C. (Bay \ BY (a) ae (34) 


where d_ is a correction coefficient which is a function of sectional area and beam- 
draft Bic derived by Vossers i960)’. The values of ie calculated by equation 
(34) are compared with those determined by Vusts! iielapleeze lhe larcve ditter— 
ence between the calculated and the experimental values should be noted. Muga's 
assumption of transverse damping force times the moment arm is not very valid. 


Therefore, the values determined by Vugts will be used in the solution of the equa- 


tion of motion. 


A3 Hydrostatic Coefficients 
The hydrostatic coefficients are easily calculable and are independent of 


frequency. 


A38.1 Heave, Z and Z 
Z © 


The hydrostatic heave and coupling coefficients are: 


Z,=-pe\ Brdé (35) 


Zeq=pe \ BY Ede (36) 


Due to constant beam, the value of the coupling coefficient, Z is zero. 


©’ 
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A3.2 Pitch, M andM 
Z © 


The hydrostatic pitch and coupling coefficients then are given by: 


M, =pg \ Bede ~ 
— 2 38 
My pg \ BYE dé (38) 


where M, is also equal to zero 


A3.3 Roll, As 


The roll hydrostatic coefficient is given by the relation: 
K, = -p8v |GM | (39) 
or 


K,=74 {GM | (40) 


where V is the immersed volume of the barge and A is the barge displacement. 


A3.4 Values 


Calculated values for the hydrostatic coefficients are listed in Table 3 in 


dimensional form for use in solving the equations of motions. 
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APPENDIX B 


WAVE EXCITATION FORCES AND MOMENTS 


Bl General Description 


The wave surface elevation of a unidirectional train of regular waves in 


deep water is defined as: 
_ 2TTX 
7(x, t) =acos own - wt +e) (41) 
where: 


Ds Jae Angular frequency of the wave 


a = Wave amplitude 
x = Distance along direction of wave travel 
X = Wave lengti: in direction of wave travel 


€ = Phase at origin when t = 0, (zero in this case) 


The wave induced excitation forces and moments are given by: 


Bow 
i(ut ~€5) i 
F.e i= dx (42) 
J stern 
BO a 
Nee od) = = —1 xdx (43) 
J stern 


The general form of the sectional wave induced force is given by: 


dF, 
<7 = |ps BREN) net) Cot) 


+U,(%) 7 Gt) | exp (- ACD) (44) 
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where: 


2 
NG) = Sectional damping [1bf - sec/ft ] 


Ls (8) = Sectional added mass [slugs/ft] 


The specific forces and moments can now be developed in terms of the hydrodynamic 


coefficients derived in Appendix A. 


B2 Heave Excitation 


The heave wave forces will be broken down into two components, one ap-. 
proaching the barge from head-on along the x-direction and the other from beam-on 
along the y-direction. The sectional hydrodynamic coefficient is found by dividing 
the total coefficient of Table 2 by the barge length since the sectional coefficients 
are constant over the length. The limits of all integrals will be from -L/2 to L/2 
unless otherwise noted and the limits will be deleted in future equations where the 


meaning is clear. 


The general expression for the surface elevation of the waves, equation 
(41) has been modified to account for the two components of the incoming wave as 


follows: 


n(x,y,t) =acos Caren + acne = wt ) (45) 


where yx is the angle between the barge heading and the direction of the incoming 


wave, 45° in this study. Then the general expression for the heave excitation is: 


al L/2 Zs 9 4. 
Z By core seo Cw 21H) cay ‘Ws -43 -o — cos ils , ay - wt 
Wt A LB NPN * Pr 


Z 
_W. em xX , omy :. ) 
+ W LB Sin nt Pan wt = (46) 
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Integrating this equation gives: 


Z, ; 
Z = (LZ)? 4 sin Te sin eX (- oni \ —=. - uw “Ww cos wt 
w,h=45° ~ \ on eet Nh) ey ILB O° «ULB 
Bs 
- W LB Sin wt (47) 


B3 Pitch Excitation 


The general equation for the pitch excitation is: 


fh Z 
| ES ae 20x _ ) “wo. (2x ) 
mt, = [a L \ x cos (Fe wt } dx + aw L \ x sin Jr wt } dx 


ZL. 
_ 2 WwW _ 2X _ _ a) 
aw L \ X COS Py se exp ( Sar (48) 
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Integrating this equation gives: 


he = Vr, 2) sin TL -L cos nae e (- ak) 2 oe sin wt 
W Jai 2TT J2 ny aN e nN L L 
ZN, 
oT cos wt (49) 


B4 Roll Excitation 


Only two components of the roll excitation will be considered since the 
dampiis excitation is known to be quite small. The effect of head-on waves has 
been neglected. The only inertial component that will be included is due to the 
shift of the center of buoyancy from the. z-axis A a+y direction creating a mo- 
ment arm for the buoyancy of each section as the wave passes. The hydrostatic 
component used was derived by Muga’. The general expression for beam-on 


waves is: 
Zz 271 H B° Zee 
— = Cc —_— eee + — " ———— ee ° 
K pw a exp( ‘ ) \ Yop? ax \ 12 sin( wt dx (90) 


where: Y CB is the motion of tae center of buoyancy in the y-direction and § is the 


sectional area under the wave. Integration of equation (50) gives: 


2 
K ep Y a exp (- 2h BPS sin BB cos ii 





IN 2m | 217 Won won 
BPA TL 
SB Ge 2 5 Sin wt (91) 
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